Economical-environmental-technical optimal power flow solutions using a novel self-adaptive wild geese algorithm with stochastic wind and solar power

This study introduces an enhanced self-adaptive wild goose algorithm (SAWGA) for solving economical-environmental-technical optimal power flow (OPF) problems in traditional and modern energy systems. Leveraging adaptive search strategies and robust diversity capabilities, SAWGA distinguishes itself from classical WGA by incorporating four potent optimizers. The algorithm's application to optimize an OPF model on the different IEEE 30-bus and 118-bus electrical networks, featuring conventional thermal power units alongside solar photovoltaic (PV) and wind power (WT) units, addresses the rising uncertainties in operating conditions, particularly with the integration of renewable energy sources (RESs). The inherent complexity of OPF problems in electrical networks, exacerbated by the inclusion of RESs like PV and WT units, poses significant challenges. Traditional optimization algorithms struggle due to the problem's high complexity, susceptibility to local optima, and numerous continuous and discrete decision parameters. The study's simulation results underscore the efficacy of SAWGA in achieving optimal solutions for OPF, notably reducing overall fuel consumption costs in a faster and more efficient convergence. Noteworthy attributes of SAWGA include its remarkable capabilities in optimizing various objective functions, effective management of OPF challenges, and consistent outperformance compared to traditional WGA and other modern algorithms. The method exhibits a robust ability to achieve global or nearly global optimal settings for decision parameters, emphasizing its superiority in total cost reduction and rapid convergence.


Motivation and incitement
The safe, secure, and inexpensive functioning of the electrical network is referred to as the optimum power flow (OPF), and this is made possible by properly configuring the system's control variables.According to the mathematical formulation, OPF is a static, nonlinear, large-scale, restricted issue with many discrete and continuous decision parameters.OPF is a vital instrument for electrical networks' cheaper and safer functioning 1,2 .OPF is used in the power system to achieve various goals, including lowering overall generating costs of thermal units (Ths), improving voltage stability, minimizing network losses, cutting CO2 emissions, and preserving the best settings for control variables 3,4 .This method is often restricted by a number of requirements that must be met, including the different technical limitations, transmission line capacity, bus voltage, ability of the power generator, and several physical limitations 5 .

Literature review
In recent literature, numerous methods of incorporating RESs into OPF and different solution methodologies for the OPF issue have been proposed.Most scholars explored bio-inspired optimization met heuristics for OPF settings in recent studies.To tackle OPF in the electrical networks with RESs, Ullah et al. developed a combination of gravitational search algorithm (GSA) and an improved particle swarm optimization (PSO) 16 , the coronavirus herd immunity optimizer (CHIO) 17 for the technical-environmental-economical dispatch problem in the two standards IEEE electrical networks.The five goals of this investigation were total fuel expenses, active losses, pollution level, stability, and variation of the network's voltage.Using the analytical hierarchy process and a weighted collect method, the multi-objective OPF was reduced to a normalized one-objective OPF.The outcomes showed that CHIO outperforms the other EETD problem-solving methodologies.Elattar concentrated on applying modified moth swarm optimization (MMSO) to mathematically model the OPF with a hybrid power and heat system through stochastic WT, and the resulting methodology was applied to the IEEE 30-bus electrical network with different test settings.In contrast to previous algorithms, RESs into OPF and the proposed technique to address it generated successful optimal solutions 18 .The OPF issue was modeled by a modified colliding bodies optimization (ICBO) 19 , and applications were made using 16 case studies on three IEEE standard electrical networks to evaluate the effectiveness and resilience of ICBO.A new, improved chaotic PSO for solving MOOPF in a test system with RESs suggested in 20 , that this PSO, compared to recent algorithms, generated superior optimal solutions 20 .In 21 , the dynamic OPF (DOPF) has been optimized by a new enhanced honey bee mating algorithm while taking into account the valve-point effects (VPEs) in the 14, 30, and 118-bus standard electrical networks.Salkuti used a new glowworm swarm algorithm (GSA) that was successfully executed to resolve an MOOPF in the IEEE systems incorporating a WT under various operating circumstances 22 .A more effective improved manta ray foraging algorithm (IMOMRFO) 23 , twenty-four benchmark tasks of varying sorts and degrees of difficulty were used to examine the effectiveness of the created approach and compare it to competing methods.The analysis's findings demonstrated that IMOMRFO produced competitive results on several MOOPFs and identified the optimum results in the recent studies for MOOPFs that arise in practice 23 .Using thermal with RESs, Kathiravan et al. studied OPF using a flower pollination algorithm (FPA).In various test situations, FPA has been examined on Indian utility 30-bus and IEEE 30-bus electrical networks 24 .The suggested solution strategy for the OPF issues has been applied to the new IEEE 30-bus electrical network using a hybrid PSO-GWO algorithm, which combines PSO with grey wolf optimization (GWO) 25 .The hybrid PSO-GWO method performed well in comparison to other algorithms, according to simulation findings, and showed that it may be a good option for solving OPF issues 25 .With adjustable WT and PV energy systems, Duman et al. solved the OPF issue using the differential evolutionary (DE) PSO (DEPSO).For testing the method under various goal functions, IEEE 30-, 57-, and 118-bus electrical networks have been studied 26 .In 27 , a new hybrid firefly-bat algorithm (HFBA-COFS) directed to handle the strictly-constrained MOOPF, which the proposed method increased the original system's capacity for population diversity and global exploration of the three different IEEE test systems to demonstrate the significant benefits of the HFBA-COFS method, the HFBA-COFS algorithm was capable of producing highquality optimal solutions, which is crucial for realizing the secure and efficient operation of massive electrical networks.In a different research, they suggested using FACTS devices to address the OPF issue.They used a modified chaotic PSOGSA to account for the uncertainty of wind energy systems 28 .A unique adaptive Gaussian teaching-learning-based optimization (TLBO) (AGTLBO) that enhanced the performance of traditional TLBO and addressed the OPF issue was proposed in 29 .The results demonstrated that it met the heuristics described in the recent studies compared to modern optimization.The AGTLBO was more efficient and successful.TLBO augmented with Lévy mutation (LTLBO) 30 was examined, assessed, and compared to other approaches using the IEEE 30-bus and IEEE 57-bus electrical networks with various OPF functions.A modified population external optimization method (CMOPEO) was suggested by Chen et al. as an improved optimizer to OPF with RESs, and CMOPEO was evaluated on the IEEE 30-bus electrical network for several test scenarios 31 .The goal function of the slime mold algorithm (SMA) 32 was the system's total cost, which included a penalty cost for underestimating RESs and a reserve cost for overestimation.Algeria's DZA 114-bus and IEEE 30-bus electrical networks were used to assess SMA performance.Four optimization strategies were compared to the SMA.According to the overall simulation findings, SMA outperformed the other analyzed algorithms over a variety of function landscapes.The 5-bus and new IEEE 30-bus electrical networks with and without unified power flow controller (UPFC) were used to evaluate the performance of bat optimization algorithm (BA)-based OPF as an effective and robust solution 33 .The performance of the turbulent flow of a water-based optimizer (TFWO) 34 to OPF was demonstrated by a comparison of the statistical indices, convergence trend, and optimal solutions to modern optimizers in recent studies 35,36 .The study's findings led the authors to conclude that TFWO was better and more successful in solving OPF optimization issues.Compared to other well-known algorithms, it had higher convergence rates and made significant technological and financial advancements.The suggested TFWO for the large-scale tested system decreases the range from 4.6% to 33.12%.The proposed solution technique resulted in a more competitive solution for the evaluated system with a notable improvement in the techno-economic aspects.Three separate target functions in OPF, including reducing overall operating costs, carbon pollution with tax, and power losses, were taken into consideration by the chaotic bonobo optimizer (CBO) 37 .To demonstrate the potency and advantage of CBO to arrive at the best result, CBO was tested on two standard IEEE electrical networks.The founded optimal solutions demonstrated the effectiveness and dependability of CBO for solving OPF using stochastic RESs, a hybrid DE, symbiotic organisms search (SOS) 38 , etc.
In the realm of optimal power flow (OPF) solutions, the landscape has witnessed a surge in innovative methodologies, each striving to address the complex challenges posed by integrating renewable energy sources (RESs) into electrical networks.The advantages and disadvantages inherent in the current state-of-the-art optimization methods employed for OPF can be summarized as follows: Advantages of state-of-the-art methods in OPF 1. Versatility in Problem Solving: • Advantage: Many state-of-the-art algorithms demonstrate adaptability and effectiveness in addressing specific challenges, such as technical-environmental-economical dispatch problems, making them versatile solutions.

Improved Performance Metrics:
• Advantage: Several methodologies exhibit superior performance in optimizing key metrics, including total fuel expenses, active losses, pollution levels, network stability, and voltage variations, showcasing their effectiveness in enhancing system efficiency.

Optimal Solutions for Specific Scenarios:
• Advantage: Certain algorithms excel in solving particular OPF scenarios, such as hybrid power and heat systems with stochastic renewable sources like wind turbines (WTs), providing tailored solutions for specific applications.

Disadvantages of state-of-the-art methods in OPF
1. Limited Scope of Application: • Disadvantage: Some algorithms may have a narrow focus, limiting their applicability to specific IEEE electrical networks or certain types of OPF issues.

Reduced Versatility:
• Disadvantage: While effective for targeted challenges, certain methodologies might lack the versatility needed to address a broad spectrum of OPF problems, potentially hindering their widespread adoption.

Limited Evaluation on Broader OPF Issues:
• Disadvantage: Some algorithms may have limited evaluations beyond specific OPF scenarios, leaving uncertainties about their performance in addressing broader optimization challenges.

Contribution and paper organization
Building upon the current state-of-the-art, our work introduces a novel self-adaptive wild goose algorithm (SAWGA) to tackle the OPF problem in a modified IEEE 30-bus electrical network with stochastic RESs, including photovoltaic (PV) and WT units.With the suggested SAWGA method and a few additional optimization techniques, we attempted to solve OPF in a modified IEEE 30-bus and 118-bus electrical networks 9 , including stochastic RESs in this work.The suggested SAWGA algorithm aimed to enhance the power of basic WGA's exploitation and exploration to resolve various operational test cases of OPF, including stochastic RESs.We have accomplished this by carefully choosing the control settings and algorithm coefficients.Lognormal and Weibull probability distribution functions were applied to stochastic RES irradiation conditions of the systems, respectively.The obtained optimal solutions by SAWGA's five evaluation items, which had been used for various OPF functions, were compared to those from the original WGA, thermal exchange optimization

Formulation of OPF in the electrical network with stochastic RESs
Traditional OPF issue was identified as an essential and vital challenge to analyze, design, and manage the power grids and energy production and transmission networks that strive to achieve the affordable and lowest cost of energy production and transmission With the condition of complying with the various stipulations and different demands, and with inequality and equality limitations of the electrical energy production and transmission system and network 9,43 .Solving this problem with various optimization functions and considering stochastic RESs such as WTs and PVs through a new modified algorithm called self-adaptive wild geese algorithm (SAWGA) has been the main goal of this research.The problem's suggested solution and its broad organization are represented below.The definition of the standard OPF issue is as follows 10 . Minimalize: Subject to: when f obj (x, u) is the OPFs' optimization function, u and x show the independent and dependent decision parameters and the equality and inequality constraints have been shown through a(x, u) and b(x, u) , respectively.

Dependent parameters
The dependent parameters in OPF have been shown in Eq. (3) 10 where P Th 1 shows the slack generator's active power, V L = V L 1 , . . ., V L NPQ shows the PQ buses' voltage values, . ., Q PV NPV are the reactive power of Ths, WTs, and PVs, and S L = S L 1 , . . ., S L NTL shows the network lines' transition power; NPQ, NTHG, NW, NPV, and NTL show the numbers of the test power system's PQ buss, Ths, WTs , PVs, and network lines.

Independent decision parameters
These parameters of OPF have been shown follows 10 If V G = V G 1 , . . ., V G NG denotes the values of all generator buses' voltage including PVs, Ths, WTs, and also, P Th = P Th 2 , . . ., P Th NTHG , P PVS = P PVS 1 , . . ., P PVS NPV , and P WS = P WS 1 , . . ., P WS NW denote the active power of Th, PV, WT units excluding the slack generator, respectively.The number of Th, PV, WT units is known as NG, NPV , and NW, respectively. (1)

Modeling of the test electrical network
This part of the study models how Th, PV, and WT power generation units may be integrated into contemporary power systems.Table 1 lists the specifications of the studied IEEE 30-bus electrical network 9 .

OPF model of Ths
According to the produced output active power, Eq. ( 5) describes the traditional OPF function in Ths as a quadratic objective function.Also, Eq. ( 6) contains the OPF model for Ths that includes VPEs, where r i and p i stand for the VPEs specifications while o i , n i , and m i are the OPF function specifications for the ith Th generator 10,45 .

Pollution level model with carbon tax C tax
Mathematically, Eq. (7) 9,10 describes how to calculate the overall emission value from heat-producing units utilizing fossil fuel.Additionally, a carbon tax model was taken into account owing to the increasing threat of global warming, and the emission cost value is determined through the addition C tax to the overall emission value that has been given in Eq. (8).In this case, CE and FE stand in for the emission costs, and sum emissions, respectively 9,10 .
where σ i , τ i , ω i , β i , and µ i are the emission characteristics of the ith Th generator.

Prohibited operating zones (POZs)
The following characteristics of a Th generator that uses fossil fuel and has several POZs 38 : where Z shows the number of POZs, z i shows the sum of POZs, P L Th i ,Z , and P U Th i ,Z−1 indicate the lower and upper bounds of the (Z − 1) th POZ of the ith unit.Figure 1 shows the specification curves of the OPF objective func- tion without and with VPEs (a), and with the POZs (b) of the Th generators.

Direct objective function models of RESs include PVs and WTs
A linear function of planned power may be used to represent a wind energy source's direct cost model 9,10 .Where DC WP i , wsh i , and P WS i stand for the ith wind power system's planned power, the direct objective function of WT, and the direct objective function specification.
(5) Where DC PV i , pvsh i , and P PVSi have been recognized as the direct objective function of PV, the direct objective function specification, and the planned output power of PV i 9,10 .

Uncertainty objective function models of PVs and WTs
The uncertain objective function models of PVs and WTs are characterized as the little-known and underestimation scenarios of these studied energy sources.The uncertainty objective function models of WTs are described in Eqs. ( 12) and ( 13) 9,10 .where UC WP,i and OC WP,i are the underestimation and overestimation objective function values, C Ow,i and C Uw,i have been defined as the uncertainty objective function specification.Also, P WS i and P wr,i are the available and rated power of WT i .A proper technique in Refs. 9,10had been applied to achieve the comprehensive formulation of the uncertainty objective functions of a PV, with the over-and under-estimation specifications C WP,i and UC WP,i , have been given in Eq. ( 3): where C Opv,i and C Upv,i represent the uncertainty objective function specification, and P PVav,i is the available or output power in the i th PV (i.e., PVi).

OPF optimization functions
Fuel cost optimization function considering VPEs Equation ( 16) models the OPF optimization function that includes the sum of the fuel costs for Ths considering VPEs and the cost of operation and transmission of the power of the WT and PV units 9,10 .( 10) This OPF problem has been modeled in Eq. ( 17).

Fuel cost optimization function considering POZs
Equation (18), which is the main goal of the mathematical model of this OPF, has been suggested to be the conventional OPF objective function in Th generators considering POZs 9,10 .

Network losses
The objective function presented in Eq. ( 19) may be used to explain the minimization of the network power losses 9,10 .
where δ ij shows the difference between voltage angles in buses i and j, and G s(i,j) shows conductance in the sth network line between buses i and j .Also, the voltage value at the ith network bus is considered by the letter V i .

Voltage deviation (V.D.)
In the considered OPF objective function, the V.D. value for the developed energy networks is determined as indicated in Eq. ( 20) 9,10 .

Equality OPF's constraints
The suggested OPF problem's equality restrictions may be formulated as given in Eq. ( 21) and Eq. ( 22) 9,10 : where P Gi shows the output power of the ith unit, P Di shows the demanded power of the ith load bus, the Q Gi shows the output reactive power of the ith generator, and Q SHi shows the output reactive power of the ith paral- lel reactive power compensator.Q Di is the demanded reactive power of the ith load bus and N bus indicates the number of network buses.The transmission line's conductance and susceptance values are denoted by the letters G ij and B ij .

Inequality OPF's constraints
• Generator limits As indicated in Eq. ( 23), there are down and up restrictions on the output power levels as well as the voltage magnitudes of the Ths, WTs, and PVs 9,10 .
• Security constraints ( 17) www.nature.com/scientificreports/Each PQ bus' voltage magnitude value must fall within certain bounds, and each transmission line's apparent power value may be constrained by its maximum capacity.Where S L i ,max and S L i indicate the maximum and apparent power levels in the i th line; V L i ,min and V L i ,max represent the minimum and maximum voltage levels of the ith PQ bus 9,10 .
where VPQ , Pslack , QTH , QWS , QPV , and SL represent the penalty factors 9,10 and "limit values'' V lim L i , P lim Th 1 , Q lim Th i , Q lim WS i , Q lim PV i , and S lim L i are defined in the same way as described in the following identity

Power and uncertainty models in the PVs and WTs
Weibull probability density function (PDF), as stated in Eq. ( 26), identifies the wind speed distribution, where ψ, , ξ , and v w indicate the scale factor, and the shape factor, the wind speed, respectively 9,10 .
Weibull fittings' outcomes in wind frequency distributions are shown in Fig. 2 9 .They were derived using a Monte Carlo simulation with 8000 iterations 9,10 .WTs' output power has been modeled as in Eq. ( 27): where the rated power, cut-in, cut-out, and rated wind speeds are denoted by p wr , v w,in , v w,out , and v w,r , respec- tively, according to wind speeds, a wind farm's electricity is divided into distinct portions, as shown in Eq. (27).The probability values are shown between Eqs. ( 28) to (30) in these sections 9,10 .www.nature.com/scientificreports/Table 2 shows the PDF components of the PVs and WTs.Wind speeds for units had been selected as v w,out = 25 m/s, v w,r = 16 m/s, and v w,in = 3 m/s, respectively 9,10 , and with selected rated power equal to 3 MW.
The Lognormal PDF was used to characterize the generation power of PV like a mathematical model of solar irradiation.According to Eq. ( 31) and Eq.(32) 9,10 , it is possible to identify the generation power and likely function of PV mathematically.
where and ξ show the standard deviation and mean of the Lognormal PDF 9 , as shown in Table 2.
By an 8000-generation Monte Carlo simulation, Fig. 3 gives the frequency distribution and lognormal probability for solar irradiation 9 .
where G pvstd , G pv and P PVrate are the standard solar irradiance level, the solar irradiance's probability level, and the rated power of PV, respectively, that have been chosen equal to 800 W/m 2 and 50 MW at bus 13.R C has been set equal to 120 W/m 2 .The generation power of PV has been shown via the histogram in Fig. 4 9 .The line shows the anticipated amount of electricity power that this PV will provide to the electrical network.It's important to keep in mind that the value of the anticipated production of solar electricity might alter.

The proposed algorithm
In this chapter, we will focus on the construction of a new SAWGA optimizer, which is based on the improvement of the previously proposed WGA algorithm.

Wild Goose algorithm (WGA)
In 49 , an effective method for high-dimension complicated and difficult optimization problems was introduced.Its name is the WGA, and it is inspired by the lives of wild geese, such as their altogether regular movement, (29) Table 2. Specifications of the WT and PV units.breeding, and development, as well as mortality in a group of geese.In the rest of this subsection, we recap the steps of the WGA algorithm.

An altogether regular migration
In PSO, each wild goose has been first identified as a population member; hence, the i th wild goose is equivalent to the i th member with the changing location X i .The migration velocity X i and individual optimal solution P i are calculated.Then, all Wild Goose populations are ranked according to their intended function, from best to worst.The wild goose migration is a collective, organized, and controlled migration that relies on power and places the upfront members so that their locations will be regularized from the most optimal solution to the worst solution (from 1 to Np ), as f P Iter i+2 ≤ f P Iter i+1 ≤ f P Iter i ≤ f P Iter i−1 .Equations ( 33) and ( 34) describe velocity and displacement in accordance with the altogether regular motion of the geese.According to Eq. ( 33), each wild goose's velocity and position changes are dependent on the speed of their front and rear goose, or V Iter i+1 − V Iter i−1 , as well as on their locations, which have been decided through the value of the optimization function of the ordered group members f X Iter i .In other words, the forefront member is the most efficient member of the wild goose since they benchmark and navigate their group as a whole as well as coordinate with the whole group.So, the migration velocity equation of the ith member by 49 has the following form: Birds in a flock always travel in the direction of the leader.The leader is in the best position globally, thus if the leader deviates, then all other birds will deviate as well.As a result, on the one hand, the best wild goose G serves as the group's leader-member, directing all other wild goose to their objective, and in other words, the position changes of these wild goose rely on their velocity.
Therefore, the amounts of displacement and change X i resulting from wild goose migration, i.e., X V i , are calculated based on the member's best local position P i and migration velocity V i , as well as the member who is in front of them P i+1 , the group leader G which has been explained in Eq. (34).The geese constantly adjust their separation from one another to maintain a safe spacing.The product of two random numbers between 0 and 1,r 7,d • r 8,d has been utilized since the quantity is small 49 :

Search for food
The ith wild goose (or any other goose) moves toward its upfront wild goose P Iter i+1 , and the upfront goose's walking and seeking for food has been modelled as its scheme.The i th member imitates the (i + 1) th member and attempts to achieve that P Iter i+1 − X Iter i .The wild geese's X W i equation for foraging while on the move is as follows 49 Reproduction and evolution of wild gooses One part of wild gooses' living similar to other alive extant has been based on the reproduction and evolution.Its modeling has been applied similarly to that of DE algorithm (Eq.( 36)), where a combination (crossover www.nature.com/scientificreports/operation) between migration equation X V i , search for food X W i has been applied.In 49 , Cr size in the WGA algorithm has been considered equal to 0.5.

Regular evolution, migration and death
The algorithm starts out in this phase with the maximum population Np initial , and as it iterates, the weaker individuals in terms of Eq. ( 37) will be eliminated from the total members (death of the weaker individuals) and the act of optimization will continue until the final members in the final generation achieves its final size Np final .Algorithm 1 illustrates the code for this simple and effective approach 49 .
where FEs max shows the maximum fitness evaluations and FEs shows the recent fitness evaluations.Sort all individuals from the most optimal solution to the worst solution.
In addition, the optimization process of the proposed modified algorithm has been shown in Fig. 5. Choosing appropriate control parameter values is typically a challenge that depends on the type of issue.Multiple optimization runs are needed since the control settings are tuned by trial and error 50 .We provide a selfadaptive technique for control settings in this section.Values for the parameter are extended to each member of the population.For Eqs. (34) to (36), we chose control-based random coefficients and modeled them for this investigation.We have recreated the equations once again and identified an appropriate connection for them, as indicated in Eqs.(38) to (40).
Cr and F are the control variables that will be modified through evolution.Both are used on an individual basis.Superior individuals result from these (encoded) control parameter values because they are more likely to live, procreate, and spread these superior parameter values to new individuals.

SAWGA for solving the different OPF problems with stochastic wind and solar power
On the IEEE 30-bus electrical network, as given in Fig. 6, the WGA, GOA, TEO, HHO, HBA, and SAWGA approaches were investigated to resolve an OPF issue, including PV and VT units.Additionally, Table 3 provides the control settings for the optimization techniques that were collected from references.The system parameters for the IEEE 30-bus electrical network have been shown in Table 4 and were derived from references 3,9,10,51 .Table 5 displays the underestimation, overestimation, and direct cost factors for PVs and WTs.The load flow framework shown in MATPOWER 51,52 has been applied to this OPF, including Ths, PVs, and WTs.For the optimization functions of the suggested OPF, all optimization methods were performed 30 times to statistically  www.nature.com/scientificreports/analyze the obtained optimal solutions.The optimization functions of the suggested OPF have been optimized in accordance with the different scenarios that follow.
• Case 1: OPF considering total cost in Th units considering VPEs, PVs, and WTs.
• Case 2: OPF considering total cost in Th units considering emission and carbon tax, and PVs and WTs.
• Case 3: OPF considering total cost in Th units considering POZs, PVs, and WTs.
• Case 4: OPF considering power losses of the electrical network.
• Case 5: OPF considering the voltage deviation (V.D.) of the electrical network.

Case 1: Total cost considering VPEs of Th units
OPF's objective function, considering VPEs of Th units, as well as the cost function of the PVs and WTs, are used in Case 1 to optimize for lowering the overall cost.In this instance, produced electric power from all various kinds of generation units employed in the redesigned electrical network are optimized to reduce the basic power cost to its absolute lowest.Table 6 lists the outcomes of all factors that were ideally determined, including reactive powers, total cost, decision parameters, and other variables-values the best, worst, average, and Std. for the algorithms and the suggested SAWGA approach after 30 separate runs are shown in Table 6.The simulation findings of this case demonstrate the efficacy of SAWGA.SAWGA offers a quick convergence trend and improved results quality compared to the traditional WGA, GOA, TEO, HHO, and HBA optimization approaches for OPF.SAWGA obtained the lowest overall power cost of 782.0238 $/h among all used algorithms.The convergences of used OPF optimization approaches have been given in Fig. 7.

Case 2: Total cost considering emission and carbon tax in Th units
By imposing a set C tax on Th units due to their CO2 emissions, the case's goal is to reduce the electricity's total cost 9 .The mandated C tax is set at $20 per ton 9 .The simulation findings clearly support the idea that imposing the carbon fee would increase the amount of wind and solar energy that is incorporated into the electricity system.

Pareto optimization method
In the realm of multi-objective problems (MOPs), the optimization process encompasses the simultaneous optimization of multiple independent objective functions (OFs), all while adhering to a spectrum of equality and inequality constraints.The formulation of these problems is as follows 39,40 : The Pareto method endeavors to identify a collection of solutions that takes into account all OFs and strikes a balance among them.The ultimate outcome of this process is identifying of Pareto optimal points-solutions that are not surpassed by any other solutions.In a broader sense, solution X 1 is considered to dominate solution X 2 if the following conditions are satisfied: www.nature.com/scientificreports/ The group of solutions not dominated by other solutions is termed "dominant solutions." To manage the size of this set, the paper employs the fuzzy grouping method.

Fuzzy grouping method
In certain MOPs, there is a need to standardize OFs to bring their values into a comparable range.Subsequently, the normalized values of OFs are merged.The fuzzy grouping method introduces a membership function (MF) for each OF, characterized by: where J min i and J max i show the minimum and maximum values of the i th OF, respectively, these values are cal- culated by optimizing each OF as a single objective problem.
For available solutions in the stored set, the normalized value of MF is computed using (44): where q shows the number of non-dominated solutions and φ h is the weighing factor pertaining to the h-th OF.

MOSAWGA to solve multi-objective OPF problems
To solve multi-objective OPF problems using MOSAWGA, the subsequent procedure is followed:     Step 1: Enter the necessary data for the algorithm and system Step 2: Embedding constraints into objective functions using the penalty functions method

Numerical results and comparison
To evaluate the performance and efficacy of MOSAWGA, it is examined on the IEEE 30-bus power system, demonstrated in Fig. 12 1,2,6 .This system has six generating units located at buses 1, 2, 5, 8, 11, and 13 and four tap-adjusting transformer units at branches 6-9, 6-10, 4-12, and 28-27.The total system power demand is 2.834 p.u. with base power equal to 100 MVA base.The lower limits of voltage magnitudes of generator buses, voltage magnitudes of load buses, VAR compensation, and transformers' tap settings are set to 0.95, 0.95, 0.0, and 0.95 p.u., respectively.Also, their higher limits are set to 1.1, 1.05, 0.3, and 1.1 p.u., respectively.
Setting the parameters of the algorithms is the same as in the fifth section.To demonstrate the efficiency of the proposed algorithm, the following three cases were taken into account: Case 1: Minimization of fuel cost and real power transmission losses.

Case 3: Minimization of fuel cost and emissions
The best solutions found by MOSAWGA in 50 trials for Case 3 are shown in Table 12.Also, Table 14 yields a comparison of all algorithms.It is deduced by investigating these tables that MOSAWGA found a solution with fuel cost equal to 862.923 $/h and emission equal to 0.2244 ton/h, which are less than those of other algorithms.The non-dominated solutions of this case are given in Fig. 15.
Case 1: The quadratic cost function for conventional generators in OPF, excluding solar and wind energy sources, can be expressed as follows: In Table 15, the performance of the proposed SAWGA is juxtaposed with outcomes from alternative algorithms explored in the field of electrical engineering.A comprehensive literature review also includes various techniques utilized for solving large-scale OPF problems.The comparative analysis in these tables showcases the superiority of SAWGA over other optimization methodologies in achieving optimal OPF solutions.The simulation results reveal a noteworthy reduction in cost, with SAWGA achieving a minimum cost of $129,522.5891per hour, surpassing the outcomes produced by alternative algorithms.
Case 2: OPF incorporating a quadratic cost function for conventional generators, along with the integration of solar and wind energy sources, can be articulated as follows: Addressing the OPF challenge involves formulating a quadratic cost function for traditional generators, accounting for their operational costs.Additionally, the inclusion of solar and wind energy sources in this scenario introduces complexities related to their intermittent nature and variable outputs.The overarching goal is to optimize the power flow in the system while considering the unique characteristics and cost implications associated with both traditional and renewable energy sources.This system is similar to the previous case study, incorporating renewable energy sources at various buses.Wind energy sources are placed at buses 18, 32, 36, 55, 104, and 110, while solar energy generation units are located at buses 6, 15, and 34.
The optimal solution for this case, obtained through the proposed SAWGA algorithm, is presented in Table 16 with a comparative study between the results of the algorithms WGA and the solutions obtained in reference 26 .These results demonstrate that SAWGA is a highly capable algorithm for optimizing and efficiently distributing loads in large, realistic power systems.www.nature.com/scientificreports/ of OPF problems in the presence of RESs, reinforcing its role as a potent tool for the sustainable and efficient operation of power systems.On the one hand, considering the simulation times in the tables, especially based on Table 15, which represents simulation results for a large-scale energy system, it is evident that the proposed method exhibits a suitable optimization speed and time, particularly compared to the original WGA algorithm.The proposed SAWGA method significantly improves the final optimal results without introducing any specific complexity or additional simulation time compared to the original WGA algorithm.

Conclusions
In conclusion, this study introduces the self-adaptive wild geese algorithm (SAWGA) as a novel and practical approach for addressing various optimal power flow (OPF) problems.We conducted a comprehensive comparison with four modern optimization algorithms (GOA, TEO, HHO, and HBA) and the traditional WGA, as well as the optimal results in the recent papers, showcasing SAWGA's superior efficiency.Our evaluation considered different objective functions within two different IEEE 30-bus and IEEE 118-bus electrical networks, incorporating wind turbine units (WTs) and solar photovoltaic units (PVs).SAWGA demonstrates remarkable capabilities in optimizing diverse objective functions and effectively managing OPF challenges.The algorithm consistently outperforms traditional WGA and other modern algorithms, achieving global or nearly global optimal settings for decision parameters.The comparison of optimization results emphasizes SAWGA's superiority in total cost reduction and fast convergence.
Our future work will focus on developing a stochastic multi-objective model for the OPF problem in the presence of renewable power generations.This avenue promises to further enhance the applicability and robustness of optimization algorithms in addressing evolving challenges in power system management.

Figure 3 .
Figure 3.The Lognormal PDF solar irradiance distribution of PV 9 .

1 :
Select the decision variables of the algorithm and the problem.2: Create the initial swarm and = 1 = [0,0, … ,0] ( = 1, 2, … , ). 3: Calculate the objective function of any member.4: Set the personal optimal location of all individuals in the group and the global best solution .5: while ≤ do 6:

Figure 10 .
Figure 10.Convergences of the optimization techniques for Case 4.

Figure 11 .
Figure 11.Convergences of the optimization techniques for Case 5.

Table 1 .
The specifications of the test electrical network.

Start Select the values of the control parameters for SAWGA. Stop condition satisfied Done Yes No Initialize a random group of geese. Evaluate the fitness the initial wild geese. Sort all individuals from the most optimal solution to the worst solution. Apply each of Eqs. (33) and (36) to each population of migratory wild geese. Is there a new member in the group which has evaluate the fitness better than that of the best member or head of the current flock of wild geese? Yes Swap the positions of that the previous best member and new member. Evaluate the fitness of the of wild gooses. Sort all individuals from the most optimal solution to the worst solution. Apply phase regular evolution, migration and death of wild gooses using equation 37. No Figure 5.
Optimization process of the proposed modified algorithm.

Table 3 .
The control variables of the optimization techniques.

Table 4 .
The pollution level and cost factors of Th units in the IEEE 30-bus electrical network.

Table 5 .
The cost factors of the PVs and WTs of the electrical network.

Table 6 .
The obtained optimal decision parameters for Case 1. Significant values are in bold.Convergences of the optimization techniques for Case 2.

Table 7 .
The obtained optimal decision parameters for Case 2. Significant values are in bold.

Table 8 .
The obtained optimal decision parameters for Case 3. Significant values are in bold.

Table 9 .
The obtained optimal decision parameters for Case 4. Significant values are in bold.

Table 10 .
The obtained optimal decision parameters for Case 5. Significant values are in bold.

Table 11 .
The best results of Case 1 by different algorithms.

Table 12 .
Best values for decision variables and the minimum fuel costs for different cases using MOSAWGA.

Case 2 :
Minimization of fuel cost and voltage magnitudes' deviance (voltage profile improvement).Case 1: Minimization of fuel cost and real power lossesThe best fuel costs and power losses of the best compromise solution (BCS) achieved by MOSAWGA and those of other algorithms are presented in Table11.The best fuel costs and real power losses achieved by MOSAWGA are 822.0989$/hand5.55MW, respectively.The final non-dominated solutions of different algorithms are demonstrated in Fig.13(only 20 points are presented).

Table 13 .
The best results of Case 2 by different algorithms.Significant values are in bold.

Table 14 .
The best results of Case 3 by different algorithms.Significant values are in bold.Case 2: Minimization of fuel cost and voltage magnitudes deviationTable12gives the summarized optimization results obtained by MOSAWGA, and Table13compares the best result of MOSAWGA with that of other metaheuristic methods mentioned in the introductory part of the article.The Pareto optimal solutions acquired by WGA-based methods are presented in Fig.14.It is observed in the figure that MOSAWGA gives well-distributed solutions compared with other algorithms.

Table 15 .
Final non-dominated solutions found for Case 3. Optimal results for Case 1. Significant values are in bold.